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Abstract 

Statistical properties of Wigner delay times and the effect of evanescent 
modes on the deterministic scattering of an electron matter wave from a 
classically chaotic 2-d electron waveguide are studied for the case of 2, 
6, and 16 propagating modes. Deterministic reaction matrix theory for 
this system is generalized to include the effect of evanescent modes on the 
scattering process. The statistical properties of the Wigner delay times 
for the deterministic scattering process are compared to the predictions 
of random reaction matrix theory. 

PACS numbers : 05.45.Mt, 05.60.Gg, 73.23.Ad, 73.50.Bk 

1 Introduction 

In the 1950's it was observed that nuclear scattering processes can have statisti- 
cal properties indistinguishable from random scattering processes jl). The first 
hint that these random elements in the nuclear scattering data might be due to 
underlying chaos in the nuclear dynamics appeared in a paper by McDonald and 
Kauffman Q who studied the energy level statistics for closed quantum billiards 
whose classical counterparts are either integrable or chaotic. They found that 
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the quantized energy levels of the chaotic billiard had a statistical distribution 
which matched predictions of random matrix theory. The first studies of the 
scattering properties of completely chaotic quantum systems with few degrees 
of freedom were due to Smilansky || and since then a number of papers have 
appeared analysing quantum scattering using semi-classical techniques 

[§,§,§, and focused on the semiclassical regime. Recently Akguc and Reichl 
studied deterministic quantum scattering from a chaotic billiard, in a regime 
where only a few channels are open, using finite elements techniques and have 
found random signatures in the Wigner delay times. 

The analysis of fully quantum mechanical scattering processes, in systems 
where only a few channels are open, is not easily accessible because this regime 
is numerically demanding. This fact has lead to renewed interest in the reac- 
tion matrix formulation of scattering theory that was developed by Wigner and 
Eisenbud [[To] in the late 1940's The idea behind reaction matrix theory is to 
decompose configuration space into a reaction region (cavity) and an asymptotic 
scattering region (lead). The exact wavefunction in the reaction region can be 
expanded in terms of any convenient complete set of states with fixed boundary 
conditions on the surface of the reaction region, provided the coupling between 
the reaction region (cavity) and asymptotic scattering region is singular 
fll3| . Reaction matrix theory provides a convenient framework for predicting 
the scattering properties of systems governed by random Hamiltonian matrices. 
We shall call the theory that uses reaction matrix theory to predict the scatter- 
ing properties of systems with Gaussian random Hamiltonians, random reaction 
matrix theory or RRMT. The predictions of RRMT have been compared to ex- 
perimental nuclear scattering data fl5f| , scattering in electron waveguides fllpf , 
and resonances in acoustic and micro- wave resonators |T^ 1 , under conditions in 
which these systems are thought to have classically chaotic dynamics. These 
predictions, in turn, can be compared to the scattering properties of chaotic 
systems. RRMT, as it is currently formulated, neglects some possibly impor- 
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tant effects in the scattering process, namely the effect of evanescent modes and 
some of the energy dependence of resonance poles. 

In this paper, we will study the deterministic scattering of an electron in 
a two dimensional electron waveguide, which has a classically chaotic cavity 
formed by a ripple billiard connected to a lead at one end (see Figure 1). This 
type of cavity is particularly well suited to the use of reaction matrix theory, be- 
cause a simple coordinate transformation allows us to construct the basis states 
inside the cavity by diagonalizing a Hamiltonian matrix. We will generalize 
the reaction matrix theory for two dimensional waveguides to include the effect 
of evanescent modes. The effect of evanescent modes on scattering processes 
has been studied for nuclear scattering processes Q , [^9| , Q , |2l| using approx- 
imate theories. For electron waveguides, we can include these effects exactly. 
We will show that for the waveguide we consider, evanescent modes dominate 
the scattering properties of the waveguide in energy regions where new propa- 
gating channels open. We will also compare the statistical properties of Wigncr 
delay times for the deterministic waveguide scattering process to the predictions 
of RRMT. 

We begin, in Section 2, by developing the reaction matrix theory of determin- 
istic scattering in our electron waveguide, starting from a configuration space 
formulation rather than the usual eigenmode formulation, and we construct the 
Hamiltonians for the cavity (reaction region) and leads (asymptotic scattering 
region) of an electron waveguide. In Section (3) (and Appendix A) we use the 
hcrmiticity of the total Hamiltonian to compute the strength of the coupling 
between the cavity and lead. In Section (4), we derive the reaction matrix. 
In Section (5), we derive the scattering matrix. In Section (6), we describe the 
method we use to obtain a complete set of basis states for a cavity with a rippled 
wall. In Section (7) we discuss the accuracy of the reaction matrix predictions 
by comparing them with a finite element calculation. In Section (8), we discuss 
the effect of evanescent modes on the scattering process, and in Section (9) we 



3 



compare the statistical properties of the Wigncr delay times for deterministic 
scattering in the waveguide with predictions of RRMT. Finally, in Section (10), 
we make some concluding remarks. 



2 Scattering Hamiltonian 

We will consider the scattering properties of an electron with mass, m, in the 
waveguide shown in Fig. (1). The electron enters from the left with energy 
E along an infinitely long straight lead which has infinitely hard walls. The 
electron wave is reflected back to the left by an infinitely hard wall located at 
x = 0. The scattering is strongly affected by the region < x < L (the cavity) 
in which the upper wall is rippled. 

The Schrodinger equation, which describes propagation of a particle wave, 
&(x, y, t), for all times, t, is given by 

B^(x,y,t) r ft 2 ( d 2 d 2 \ . .1 . 

where ft is Planck's constant. The potential, V(x,y), has the following prop- 
erties: V(x,y) = oo for (L<x < oo); V(x,0) = oo for (— oo<x<L); V(x,y = 
g(xj) = oo for (0 < x < L); and V(x,y = a + d) = oo for (— oo < x < 0); 
where g(x) = d + acos(47ra; / L) gives the contour of the ripple, d is the average 
width of the cavity, L is the length, and a is the ripple amplitude. Throughout 
this paper, we take the electron mass to be the effective mass of an electron in 
GaAs, m = 0.067m e , where m e is the free electron mass. 

We can introduce projection operators, P = J^^dx J^dy \x,y)(x,y\ and 
Q = fo^x J^dy \x, y) (x, y\ which satisfy the completeness relation Q + P = 1 
(all wavefunctions are zero for L < x). Here \x, y) is the simultaneous eigenstate 
of position operators, x and y. The projection operators, Q and P have the 
property that Q — Q 2 , P = P 2 , and QP = PQ = 0. If a state, |*) has spatial 
dependence, ^(x,y)=(x,y\^), over the interval (— oo < x < L), then the state 
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(x,y\Q\V) = &(x,y) for (0 < x < L) and the state {x,y\P\$) = ^(x,y) for 
(-00 < x < 0). 

Inside the cavity (Region I, (0 < x < L) in Fig.(l)), we define a Hamiltonian, 

1 



Hqq=Q 



2m (p 2 x +p 2 y ) + V(x,y) 



Q, (2) 



where p x and p y are momentum operators and m is the mass of the particle. The 
Hamiltonian, Hqq, is Hermitian and therefore it will have a complete, orthonor- 
mal set of eigenstates which we denote as Q\<j>j). We can write the eigenvalue 
equation in the region, < x < L, as HqqQ^j) — XjQ\(pj), where Xj is the j 
energy eigenvalue of Hqq and j = 1, 2, ...M (we will later let M— »oo). Because 
there is an infinitely hard wall at x — L, the eigenstates 4>j(x, y)=(x, y\Q\4>j) 
must be zero at x = L. We have some freedom in choosing the boundary condi- 
tion at x = 0. In this paper, we will require that the eigenstates, <pj(x,y), have 
zero slope at x = so that ~^-\ = 0. Singular coupling, between the cavity 
and the lead, will correct for the fact that the actual wavefunction does not 
have zero slope at x — 0. The completeness of the states, Q\<j>j), allows us to 
write the completeness relation, ^2jQ\4 > j)(4>j\Q = Q- Orthonormality requires 
that (<pj\Q\4>j>) = dj.j'- The part, inside the cavity, of any state, |\E r ), in the 
waveguide can be expanded in terms of this complete set of states, so that 

M 

QI*> = E<&l6l*> l&>- ( 3 ) 

i=i 

Inside the lead (Region II, (—00 < x < 0) in Fig. (1)), we define a Hamilto- 
nian 



H P p=P 



1 ',2 



P. (4) 



Its eigenvalues are continuous and have range, (0<E<oo). The eigenvector 
of Hpp, with eigenvalue, E, will be denoted P\E). The eigenvalue equation 
then reads, H PP P\E) = EP\E). Because the leads are assumed to be straight, 
the transverse parts of the energy eigenstates, in the leads, decouple from the 
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longitudinal part. Because the walls of the channels are infinitely hard, the 
energy eigenstates in the leads (for x < 0) can be written 

N 

p\e) = ^r n p|$ n ), (5) 



n=0 



where r„ = ($ n \P\E) and 



{*MP\*n) = xkA*))l^to(^) (6) 

represents the contribution to (x\E) from the n th transverse quantum state 
in the lead. Although we have summed over the first N transverse states we 
will later let N-^oo. For a particle with energy, E, the state, P|<J?„), has the 
property that 

H PP P\* n ) = PP|*„> = ^(kl + (ifA 2 )Pl*n). (7) 



where 



The state, P\$> n ) is called the n th channel. Eq. (||) gives the decomposition of 
the total energy, E, into its longitudinal and transverse parts when the electron 
is in the channel, P|$ n ). For a given energy, E, there are an infinite number 
of channels for the particle, some propagating and some evanescent. Channels 
with propagating modes occur if the longitudinal wavevector is real. Channels 
with evanescent modes occur if the longitudinal wavevector is pure imaginary. 
Evanescent modes describe localized contributions to the electron states in the 
waveguide. There are an infinite number of them and for some values of the en- 
ergy, E, they play a dominant role in determining the dynamics in the waveguide 

H). 

We couple the cavity and the lead at their interface, x = 0, via the singular 
operator, V — CS(x)p x . The coupling constant, C, can be determined by the 
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condition that the total Hamiltonian Jul be Hermitian (see Section (3)). Then 



Hqp = QVP = C- dxxl dyi / dx \x 1 ,y 1 )5(x 1 -xo)5(x Q )- — (x ,yi\ 



dxo 



(9) 



and 



H PQ =PVQ = C- dx dy / dxi \x Q ,y Q }5(x Q - xi)5(xi)- — (xi,yo|- 

IJ-oo J -co Jo dx l 

(10) 

It is useful to remember that J Q L dx 5{x) — |, J Q L dx 5(x — xq) = 1 if < xq < L, 
and J Q L dx S(x — xq) — if L < x$ or xq < 0. Note also that 
r L 

Hqq = I dx 



-Ti 2 ( d 2 d 2 



and 



Hpp — 



dx 



-fi 2 ( d 2 d 2 

? y \ x ^W + w* +v{x > v)]M 1121 



The total Hamiltonian of the system can be written 
H = -Hqq + Hpp + Hq P + H P q. 



(13) 



The waveguide energy eigenstates, \E), satisfy the eigenvalue equation H\E) = 
E\E) . The states, \E), can be decomposed into their contributions from the 
two regions of configuration space, so that 



M N 
3 = 1 n=0 



(14) 



where jj 
the form 



i\Q\E) and T n = (<& n \P\E). The eigenvalue equation then takes 



(Hqq 
H ol 



Hpq Hpq 
Hpq Hpq 



Hqp Hq P 

Hqp Hqp 

Hpp 

Hpp 



/7iQI0i>\ 



TiP|$i) 

r 2 P|$ 2 ) 



V 



/7iQI<M\ 



E 



TxP\*i) 

r 2 P|$ 2 ) 



V 



(15) 
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This yields a series of equations 

N 

HqqQI^H + ^ff QPJ p|$„)r„ = EQ\<t>j) 7j , (16) 

n=0 

for j = 1,2, ...,M and 

H PP P\<f> n )T n + Y^PoQWiiH = E P\®n)T n (17) 

j 

for n = 1,2, ...,N. The condition for Hermiticity of the Hamiltonian, ("J^i? |$ a ) 
(*f? a \H\if? p)* , allows us to determine that the value of the coupling constant, C, 
is C = |^ (see Appendix A). 

3 The Reaction Matrix 

We now have enough information to derive the reaction matrix for this system. 
Let us first multiply Eq. ( |l6| ) by (<pj\Q to obtain 

N 

(hlHQQlhhj + ^(hWQP^Tn = E{4>j\Q\4>i)li (18) 
n=l 

which reduces to 



4i^^ n ' ' dx 

n—1 



r„ = o. (19) 





If we use Eq. (|67j) for the coupling constant, C, we can rewrite Eq. ( |19[ ) and 
obtain the following expression for 7^, 



ft 2 1 v^, t / n \ d X™ 

7 ^ = SIFj5jM» (0) ¥ 

v J ' n—1 



(20) 



The continuity equation (p5|), when applied to energy eigenstates, yields 

M N , 

r„X„(0) = ^7i0i,n(O) = 51 ^ IV (21) 

j=l n'=l 

where 

n 2 ^lAWsAo) 
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= %_aJ (22) 
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is the (n, n') th matrix element of the reaction matrix. 

We must now distinguish between propagating and evanescent modes. The 
states in the leads, for propagating modes, can be written 



where 



+ 



fan 



>2mE 



(23) 



(24) 



% A \d + a/ 

If there are v propagating modes then n — 1, 2, v. Here we use a unit current 
normalization. The evanescent modes in the leads can be written 



T nXn (x) = -^e-^W, 



(25) 



where 



d + a 



2mE 



(26) 



For evanescent modes the index n = v + 1, v + 2, N where N-^oo. 

4 The Scattering Matrix 

To obtain the scattering matrix, we must first separate the propagating modes 
from the evanescent modes. This first step is accomplished as follows. Using 



Eq. ( |23|) and Eq. (25) we can write Eq. (|2l| ) in the matrix form 



a + b 

c 



K p W R pp R pe W K p \( i(b - a) 



K e 



where 



Rep -*^ee 



K e 



(27) 



K„ 




Cjv 





Ri,i 



R 



1,1/ 



/ Ri,v+i 



R 



VP 



i Rpe 



Rep — 



Ru,l ■■■ Ru,u ) \Rv,v+\ 

f Rv+i,i ■■■ Rv+i,v\ / Rv+l,U+l 

: : ' R £e = : 

\ Rn,i ■ ■ ■ Rn.u ) \ Rn,i> 



+1 



If we expand out Eq. (E7]) , we find 



a + b = iK p R pp K p (b — a) + K p R pe K e c 



Ri. 



N 



Rv,N , 
Rv+l,N 



Rn 



N 



(28) 



(29) 



c = iK e R ep K p (b - a) + K e R ee K e c. 
From Eq. @ we can write c as 

-K e R ep K p (b - a), 



(30) 



(31) 



(l e - K e R ee K e ) 

where l e is a unit matrix with the same dimensions as R ee . If we substitute Eq. 
© into Eq. (§f), we find 



a + b = iD(b — a), 



(32) 



a e -K e R ee K e ) KeRepKf 



(33) 



where 

D= K p R pp K p + K p R pe K\ 

The second term on the right in Eq. ( J33| ) contains the effect of the evanescent 
states on the propagating modes in the waveguide. The scattering matrix, S, 
relates the outgoing propagating modes to the incoming propagating modes 
through the relation, a = Sb. The scattering matrix is thus given by 



5' 



(l P -iD) 



(34) 



(l P + iD) 

where l p is a unit matrix with the same dimension as R pp . We see from Eqs. 
( |33| ) and (|34|), that the evanescent modes may play an important role in the 
scattering process. To see this effect on the resonance structure of S matrix, we 
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obtain a more explicit form as follows. First we define the coupling matrices, 
f <f>u ... <f>l p \ / (j>n ... <f>u \ 

w Np = : : \ K p and lo^e = : : K e , (35) 

\ 4>Nl ■ ■ ■ <f>Np } \ 4>N1 ■ ■ ■ 0JVe / 

where p is the number of propagating modes and e is the number of evanescent 
modes in the lead. The matrix D can be written in terms of the coupling 
matrices as, 

3 = "U rn ft ™ Np 

1 11 

+ ™In pf S-^ e ? - t i - — ^eJV^r S-^p( 36 ) 

tLl N - H in l e - to^y ElN _ H . n W Ne tiLN-tiin 

where /?i„ is a diagonal matrix formed by the eigenvalues of the Hamiltonian, 
Hqq, inside the cavity. The second part of Eq. (|3^) can be rearranged by 
expanding in a series and regrouping terms, 

D = fj ^Np 



EIn — Hi n 

"In — 

N — tl in — WjSf e W eN 



1 ,t 1 



W vN "Ft 7^— w NeW eN — = r— W Np . 

pJV El N - H in E1 N - H, - »"-- »»t 



This expression for D can then be substituted into Eq. (34) and we obtain the 
following form for the scattering matrix, 

S = -(l-2iw f N —z = i — ^WNp) (37) 

E1 N - H in - w Ne w' eN + iw Np w pN 

As can be seen from the denominator of this expression, evanescent modes 

affect the positions of resonance poles in complex energy plane because of their 

dependence on both coupling matrices, WpN and w e ^ . Contributions from term, 

WN e w\ N are not included in RRMT calculations. Also, the energy dependence 

of the coupling matrix, w p m is neglected in RRMT calculations, although this 

is known to various authors and they simply assume that the energy regions 

they consider are far from channel openings. As we see in Eq. (^), evanescent 

modes may play an important role in the scattering process, and in subsequent 

sections, we will investigate their effect on scattering of an electron from the 

ripple cavity. 
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5 Basis States for the Cavity Region 

We now describe a method to obtain the complete set of eigenstates, Q\4>j), of 
the Hamiltonian, Hqq. We will require that these states have zero slope at the 
at the cavity- lead interface [x = 0). We introduce a coordinate transformation 
which straightens the rippled wall of the cavity || , (3(| . Then we can obtain 
a Hamiltonian matrix which can be diagonalized to find the eigenvalues, Xj, 
and eigenstates, 4>j{x,y). The first step is to write the eigenvalue equation, 
HqqQ^j) = m configuration space. It takes the form 

-h 2 ( d 2 d 2 \ 

l^\d^ 2 ^dy 2+V{ - X:V ^ jfofavi^Wjfav)' ( 38 ) 
where <j>j(x,y)=(x,y\Q\cj)j). After the coordinate change, 

u = x, v = — y 4 , (39) 

a + a cos{^x) 

we obtain an eigenvalue equation in terms of the coordinates, u and v, given by 
h 2 

5^{u,v) = - — {d 2 u + h x d 2 v + h 2 d 2 uv + h 3 d v )^{u,v) = X^ji^v) (40) 
where 

1 + v 2 gl _ -2vg u -vg uu 2vgl 

m — n j "2 — 7 ft 3 — 1 9— j 

9 9 9 9 

g = g{u) = d + a- cos(^u), g u = |f , and ipj(u,v) = <j>j(x(u,v),y(u,v)). The 

boundary conditions in (u,v) space are given by d u ipi(0,v) = 0, ipi{L,v) = 0, 

i/ji{u,0) = 0, and ipi(u, 1) = 0, so that in terms of these coordinates the walls 

are straight. Note that in the (u,v) coordinate frame, the states, tpj(u,v) are 

normalized with a weighting factor, g(u), so that 



g(u) ipj(u,v)il}ji(u,v) dudv — Sjj>, (41) 

The state, ipj(u, v), can be expanded in terms of a Fourier basis, 

00 00 

") = EE Bl n cj> mn (u, v) (42) 



m— 1 n— 1 
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with 

/ x 2 i . , , ,(2m— l)iru. , . 

<pmn(u, v) = -j=g 1/z sm(ri7™)cos(- — ), (43) 

where B 3 mn are the unknown expansion coefficients. As a result of this expan- 
sion, the boundary value problem is transformed into the eigenvalue problem, 



oo oo 



£ E Hmnm'n'Bi^, = EjB^. (44) 
m—1 n— 1 

The Hamiltonian matrix elements, H mnm i n i, are given by 

— I du dv^/g sin(nirv) f H {sin(n' itv) f 1 / \fg) (45) 



L 



o Jo 



where / = cos( t 2 '"^ 1 )™ ^ f = cos( ^ 2m 2 ^ nu ), g = d+a cos(^u), and H is the 



differential operator defined in Eq. (40). Note that we cannot use integration 



by parts to get a symmetrical form, as was done in Refs. j|] and Q because 
surface terms will not drop out. 

We calculate the Hamiltonian matrix elements using Eq. (|4^). We reduce 
the double integral to a single integral after integrating in the v direction. After 
some algebra we find the following form which is suitable for numerical calcu- 
lations, 



Hmnm'n' = ~W mm ,8 rm ' + ( n7: ) 2 w mm' ^nn' /2 + (2(mr) 2 - 3)w^ nm ,S nn i 

+< m ,<W/4 + (1 - ^nn0(-l) n+n, ( f ? f tT "'^2 + 7r' n t0 ( 46 ) 

[n l — n' A y (n z — n' z ) 



where 



(i 



i / , 4/ u „5 2 - 4f u g u g - 2fg uu g + 3/g 2 , 



w 2 mm , = [ du^f,wf nm , = [ du-^f 
Jo 9 Jo 9 

4 _ f L , -V u g u g + fgl - fgu U g + Vgj t , 

Wmm> = / du ^ / 

Jo 9 
The eigenvalues and eigenvectors of H can be calculated efficiently due to the 
sinusoidal integrals. Eigenvectors of H give values for the expansion coefficients, 
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B° mn , and the eigenfunctions in u-v space can be found from these coefficients. 
The solution can then be transformed back to x-y space to obtain the basis 
states, 4>j (x, y). 

6 Accuracy of Reaction Matrix Theory Compu- 
tations 

We have computed the scattering matrix and the amplitude of the evanescent 
modes for the waveguide given by Fig. 1, using both the reaction matrix theory 
presented in the previous sections, and an independent finite element method 
(FEM) as a check on the reaction matrix results. For this comparison we used 
the following parameter values: a = 10A, d = 100A, L = 20(M, and m = 

9 fc2 2 

0.067m e where m e is the free electron mass. This gives E n = n 2m(d+a) 2 = 
(0.0622n 2 )eV r for the energy at which the n th propagating mode appears. 

We will discuss our results for the energy interval, E\ < E < 9Ei, in which 
one propagating mode [E\ < E < AE\) and two propagating modes (4Ei < 
E < 9Ei) can exist in the leads. We studied both the case in which a = 0, so 
the cavity is rectangular, and the case a — 10A in which the dynamics in the 
cavity is fully chaotic We have computed the partial Wigner delay times, 
t„ = where 6 n is the n th eigenphase of the S-matrix. 

For a rectangular cavity, a — (the upper boundaries of the rectangular 
region is shown with a dotted line in Fig.l), we can find analytic expressions for 
the S-matrix which serve as a check on the accuracy of our programs. For the 
energy regime where only one mode can propagate in the leads, the S-matrix 
(S-function in this case) is a complex number with unit magnitude and is given 
by the value of reflection coefficient, S — e l2kL . This is the phase shift of the 
wave as it enters the cavity and reflects back to the entrance. The eigenphase 
of the S-matrix is 9 = 2kL. Since the phase angle depends linearly on k, no 
resonance occurs. The Wigner delay time has no peaks. 

The reaction matrix for a rectangular cavity can be written exactly. Since 
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there is no mode coupling in the rectangular cavity, it is enough to calculate 
Rn in the energy regime where only one mode propagates. The eigenvalues of 
Hqq take the form ^ ( ^ m d P — h ^ 2 '(^l)'^ )' wnere m an d n positive integers 
representing the transverse and longitudinal degrees of freedom. We obtain 

N 1 
Rn = lim > 



N ^*> H e " (V^) 2 - ((2 • n - 1) ■ tt/(2L))2 



(47) 



n=0 

tan(L • - (tt/cZ) 2 )) _ tan(fci) 

^e-W)) = *~ ' 

where e = 2mE/h 2 . The S function is given in terms of Rn, 

q= 1 - » fc ^n _ 1 + j tan(fcL) j2fcL 

1 + iftiZn l-itan(Jbi) ' 1 ; 

so the eigenphase is given by = 2fcL which coincides with the exact result. 

We can now check the accuracy of the methods we are using. In Fig. 2, the 
sum in Eq. ( [47j ) is truncated at N = 1000 and N = 10000 and used to compute 
the reflection coefficient. The relative error is less than 10 for the number of 
terms we kept in following calculations. The phase angle curve calculated by 
FEM is in agreement to the order of 10~ 5 with the exact result. More discussion 
about the FEM method can be found in ref j^] and references cited there. 

Let us now consider the case of the ripple cavity. In Fig. 3, we show the 
Wigner delay times calculated by using Eq. (|^) and we compare them to the 
Wigner delay times obtained from the FEM calculation. In Fig. 3(a) we show 
the Wigner delay time, n, in the energy interval where there is one propagating 
mode in the leads. In Fig. 3(b), we show two partial Wigner delay times, t\ and 
T2, in the energy interval where there are 2 modes in the lead. We have kept 
up to 2500 terms in Eq. (|2^). We checked the accuracy of these eigenvalues 
by increasing the size of the Hamiltonian matrix and comparing eigenvalues 
of a matrix with 5500 eigenvalues, and a matrix with 10,000 eigenvalues. We 
found that the first 2500 eigenvalues were the same to an accuracy greater than 
10 -5 . We used the first 2500 eigenvalues and their eigenstates to construct the 
reaction matrix. 
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We also looked at the analytic continuation of the S-matrix into the lower 
complex energy plane to see the resonances explicitly. In Fig. 4 we used Eq. 
([54]) (including evanescent modes) to obtain the analytic continuation of the 
S-matrix. The position of the three peaks, shown in Fig. (4), are E\ = 1.3585 — 
i0.123, E 2 = 1.8991 -£0.2403, and E 3 = 3.1241-^0.2316. The real part of peak 
position shows the resonance energy and the imaginary part shows its lifetime. 
These are both consistent with Fig. (3. a). The full width at half maximum 
of resonances in Fig. (3. a) is in agreement with imaginary parts of S-Matrix 
poles in Fig. (4). We obtained the following numerical values for the resonance 
positions, J?j, and their widths, Ti in Fig. (3. a), (Ei — 1.3733, Ti = 0.16), 
(E 2 = 1.8917: T 2 = 0.20), and (E 3 = 3.1359, T 3 = 0.21), where T is defined as 
the full width at half maximum. In the upper half plane one gets corresponding 
zeros of the S-Matrix. 

7 The Effect of Evanescent Modes 

For the shape of waveguide cavity that we consider here, the effect of evanescent 
modes on Wigner delay times appears to be most important just before a new 
channel opens in the lead. We have studied the effect of evanescent modes 
using parameters, a = 10A, d = IOOA, L — 500A, and m = 0.067m e , and we 
use the form of the S-matrix in Eq. ( |34|) to obtain our results. We compare 
the variation of S-matrix elements, Sij, for an S-matrix which includes the 
evanescent modes (K e j^0), with an S-matrix, S , which excludes evanescent 
modes (K e = 0). We use the cavity length, L = 500A, (rather than the samller 
length L — 200A used to check accuracy) to increase the density of resonances 
in any given energy interval. In a mesh based numerical method (like FEM or 
a finite difference) increasing cavity length is numerically is not efficient due 
to the increasing number of nodal points, but the reaction matrix approach 
can easily accomodate longer cavitites. In Fig. (5. a) we show the effect of 
the evanescent modes on the Wigner delay time at energies just below where 
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the second channel opens and, as we expect, there is a considerable increase 
in the delay of the electron. The absolute value of the amplitude, ci, of the 
first evanescent mode also increases just before the second propagating channel 
opens as shown in Fig. (5.b). The effect of the first evanescent mode, ci, is 
dominant since the amplitudes of the second and higher evanescent modes are 
near zero. A similar behavior of the evanescent modes occurs at energies just 
below where the third propagating channel opens. There again, one evanescent 
mode becomes dominant. In Figs. (5.c) and Fig. (5.d) we compare the behavior 
of first and second partial Wigner delay times, n and T2, respectively, both for 
the case when the contributions of the evanescent modes are included and for 
the case when they are removed in calculations of the S-matrix in this energy 
regime. 

In Fig. 6 we plot |5n| and 15*121 in the energy interval, 4Ei < E < 9Ei (two 
propagating modes). We find that 15*221 = |5n| and |52i| = |5i2|. Therefore we 
show only these two matrix elements. In Fig. 7, we show the effect of evanescent 
modes on S-matrix elements in the two mode energy regime by plotting the 
differences, 15^1 — |5n| and |5° 2 | — |Si2|- The difference in the magnitude of 
the S-matrix elements is small, but the difference in the slopes can be fairly 
large. 

We have also looked at the analytic continuation of S-matrix elements in 
the complex energy plane and we find good agreement with the predictions of 
Wigner delay time plots. In Fig. 8, we show partial Wigner delay times in 
energy interval, 4E\ < E < 9Ei. In Fig. 9 we show the behavior of |5n| in the 
complex energy plane. Fig. 9. a gives large scale behavior, and Fig. 9.b focuses 
on behavior near the real axis. The poles near to real energy axis (shown in Fig. 
9.b) determine the sharp peaks in the Wigner delay times. The poles further 
from the real axis determine the broader peaks in the Wigner delay time plots. 

In Fig. 10 and in Fig. 11 we show the effect of the energy dependence of the 
coupling matrices, wn p and wjy e . This energy dependence is always neglected 



17 



in RRMT calculations. In Fig. 10, we plot \Su\, both for the case when the 
energy dependence of WNp and WNe is taken into account (full line), and for the 
case when the energy dependence of wn p and WMe is fix at value, E = 6.5-Ei 
(dotted-dashed line). In Fig. 11, the effect of the energy dependence of wn p 
and wpf e on the distributions of poles in the complex energy plane is shown. 
The position of S-matrix poles changes when the variation with energy of the 
coupling constants is not included. In Fig. 11, the solid lines are contour lines 
of | Snl for the reaction matrix calculation with evanescent modes included and 
the energy dependence of wn p and WNe included as in Fig. (6). The dotted- 
dashed lines shows the same quantity but using coupling matrices, wn p and 
WMe: with dependence on energy, E fixed at the real value E = 6.5i?i. Neglect 
of the energy dependence of the coupling constants causes a shift of the poles 
away from their true positions. This shift is small in the neighborhood of the 
fixed energy, E — 6.5-Ei, but it grows as one moves further away in energy. 

8 The Signatures of Chaos 

In this section we compute the statistical properties of the Wigner delay times 
obtained for deterministic scattering of the electron from the ripple cavity. We 
consider only configurations of the ripple cavity for which the dynamics of the 
cavity is classically chaotic. One can use either Eq. (j|J) or Eq. (^) to calculate 
Wigner delay times deterministically. We have checked that they give identical 
answers. For the deterministic calculations, we can get sufficient data to develop 
good statistics by changing the ripple size from a = 1(M to a = 30A in units of 
0.2A. In Figs. (9a), (9.b), and (9.c), we show the statistics for the total Wigner 
delay times for deterministic scattering for the cases when M — 2, M — 6, and 
M = 16 propagating modes, respectively, exist in the leads. In these Figures, 
P(t) is the histogram of Wigner delay times normalized so the area is equal to 
one, and < r > is the mean Wigner delay time. The distribution, P(t), shifts 
from a Poisson-like distribution to Gaussian-like distribution as we increase 
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the number of channels. For a small number of channels the distribution is 
asymmetric and has a long tail. 

We also looked at the statistics of the total Wigner delay times obtained by 



replacing the S-matrix in Eq. (37), by the equation 

§=-(l + 2ig 2 ^ — 1 w), (49) 

.Mat — n in — ig ww< 

where H in is chosen from an Gaussian Orthogonal Ensemble (GOE) and the 
coupling matrix, w, is constructed from M eigenvectors of the realizations of 
H' in in the same GOE ensemble. Note that g is a coupling constant that must 
be determined from experiment. We also checked our result by building w using 
the M eigenvectors of each realization of Ht n and we get a similar distribution 
for the corresponding number of channels. 

We have calculated the Wigner delay time by taking the derivative of the 
S-matrix eigenphase curve, 9(E) versus E, in two different ways. The first way 
is to take two neighboring energy points (we chose E=0 and E=0.001) and used 
these obtain one Wigner delay time for each realization of H in . The second way 
is to obtain a whole seeries of Wigner delay times from the 6(E) versus E curve 
for a single realization of H in . We have checked that these two methods give 
similar results as we would expect due to ergodicity. 

In Fig. (12) we show the distribution of total Wigner delay times obtained 
from the Gaussian Orthogonal Ensemble as described above. The middle col- 
umn with Figures Figs. (12. d), (12. e), and (12. f) corresponds to a coupling 
constant, g — 1.8, which is the strong coupling regime for RRMT. The right- 
most column with Figures Figs. (12. g), (12. h), and (12. i) corresponds to a 
coupling constant g = 1.0. The distribution of total Wigner delay times for our 
deterministic scattering from the chaotic ripple cavity, agrees qualitatively with 
the predictions of random matrix theory for strong coupling. This is consistent 
with the fact that the opening between the ripple cavity and the leads for our 
case is very large. 
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9 Conclusion 



In this paper we have studied the effect that evanescent modes have on the 
scattering properties of an electron in a waveguide with a "chaotic" cavity. We 
have reformulated the reaction matrix theory of electron waveguide scattering to 
explicitly include the effect of evanescent modes. We have found that evanescent 
modes can increase the delay of the electron for energies near the opening of new 
channels. This effect has been seen before |25|]. The scattering system we have 
considered is relatively "soft" . There are no impurities and no sharp corners to 
snag evanescent modes, and yet their effect is still noticeable. For systems with 
impurities and sharp corners, we expect the effect of evanescent modes to be 
even more dramatic. 

We have also studied the effect of neglecting the energy dependence of the 
coupling matrices that appear in the reaction matrix approach to scattering. 
This appears to cause an efective repulsion on the positions of quasibound state 
poles. 

The effects of both the evanescent modes and the energy dependence of 
coupling matrices are routinely neglected in RRMT, and this should be kept in 
mind when attempting to use that theory to make predictions about real waveg- 
uide scattering experiments or numerical simulation of deterministic waveguide 
scattering systems. 

We have also studied the statistical distribution of the Wigner delay times 
for scattering from our chaotic waveguide cavity, for the case of M = 2, M = 6 
and M = 16 propagating modes. To build adequate statistics for comparison 
with RRMT predictions, we have included data for a range of ripple amplitudes, 
being careful to include data only from the regime where the internal dynamics 
of the ripple cavity is completely chaotic. If the ripple amplitude is too large or 
too small, the cavity will again develop a mixed phase space |50| . We find fairly 
good qualitative agreement with the predictions of strong coupling RRMT. 
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11 Appendix A: Hermiticity Condition 

Consider the arbitrary states, \^f a ) and l^p). The condition for Hermiticity of 
these states is that 

{*p\H\9 a )-{* a \H\9[,y = 0. (50) 

We will use this condition to determine the coupling constant C (this method 
of determining the strength of the coupling was first suggested by Pavlov |l4| ) . 
We can decompose the states \^ a ) and \^p) into their contributions to the two 
disjoint configuration space regions and write them in the form, 

M N 

|*a) = Q\*a) + P\^ a ) = &-Q|^> + 5>„P|$„), (51) 

j=l 71=1 

where aj = (<j>j\Q\i& a ) and A n = ($ n |P|* ), and 

M N 

j=l 71=1 

where bj = (<f>j\Q\typ) and B n = (& n \P\if> p) , Inside the cavity, 0<x < L, we 
have expanded \^ a ) and in terms of the complete set of energy eigenstates, 
Q\4>j), of the Hamiltonian, Hqq. In the lead, — oo<x < 0, we have expanded 
\^ a ) and l^p) in terms of the complete set of energy eigenstates, P|$ n ), of the 
Hamiltonian, Hpp. 

The Hermiticity condition takes the form 

{9(}\H\9 a ) - {* a \H\9fiy 

M 

= £[<&I#QqI&) - (4>j I Hqq I <j>j )*]dj bj 
J=l 
M N 

EE^A*) - (<S>n\Hp Q \^y}a*B n 
j'=ln=l 
M N 

+EEk $ ™i#«?i^> - (hi&QppnyfaK 

j=ln=l 

N 

+J2i(®n\Hpp\<S>n) - {$ n \Hp P \3> n )*}A* n B n = 0. (53) 
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We can now evaluate Eq. ( p3[ ) term by term. The Hamiltonian, Hqq, is 
Hermitian and its eigenvalues are real so we immediately have (4>j\HQQ\4>j) — 



\Hqq\03 



A,- — A* = 0. We will assume boundary conditions 



dx \x—a 







and use Eqs. (|^) and (J10|) . Then /or these special boundary conditions we find 
(3> n \HpQ\<f)j) = and (& n \HpQ\4>j)* — 0. Thus, the Hermiticity condition 
reduces to 



(*/j|#|*a> - (*a|^|*/j)* 



M iV 



j=ln=l 



AT 



-^[($„l^p P |$„) - ($ n |ffpp|$ n )*]A;£„ = o. 



(54) 



It is useful now to perform the spatial integrations implicit in Eq. (p4|). Let 
us first consider the last term. We can write & n (x, y)=(x, y\<& n ) — Xn( x )y' di^ sm (f+|) 
If we substitute into the last term in Eq. (|54|) and perform the integration over 
y, we obtain 



2m 



dx 



2m 



($ n \Hpp\$ n ) - ($ n \Hpp\$n)* 
d 2 d 2 
dx ax 

d d 1 

Xn(x)-1-Xn(x) ~ Xn{x)—X*n{x) 

dx dx 



(55) 



Note also that 



b } \H QP \* n ) =C- / dx x 



s(»i) 



dyi 



dxa 



dyo 



(56) 



x((f> J \xi,yi)5(x 1 )5(x 1 - x )6(yi - y ) — (x ,yo\$ n ). 

Perform the integration over x\ and notice that 17(0) = d + a. Then Eq. ( |56| 
takes the form 



t> j \H QP \3>„} = C- 



c£+a 



dyi 



dxQ 



d-\-a 



dyo 



x(<Pj\0,yi)S(x )S(yi - y ) — (x ,yo\® n ) 



(57) 
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The cavity basis states, at the interface, can be written 

nny 



(&|0,tf> = #, n (0)^ 



-sin 



d + a \d + a 

Thus if we perform the remaining integrations in Eq. (p7[), we finally obtain 



(58) 



6 J |ff O p|$ n )=C-^ jn (0). 



fir: 



(59) 



(i 



We can now combine the above results and write the Hermiticity condition in 
the form 

_ M N 



j=ln=l 



^iEE(*;...(o)^ 



dx 



a*B n - 0j,n(Q)~ZT 



2 iV 



— ft x - 

2m ^ 



A*B n = 



(60) 



for the case of zero-slope boundary conditions for the cavity basis states. The 
boundary conditions at x = — oo cannot depend on details of the interface at 
x = 0. Therefore we must satisfy the conditions, 

_ M N 



j=ln=l 

-2 N r 



fix 



2m 



E 

n=l 



and 



-2 JV 



— ft x - 

2m ^ 



X*n{x)^Xn{x) - Xn(x)-^X*n{x) 



X*n(x)-^Xn(x) ~ Xn(x)^-X*n(x) 



A* n B n = 



x=0 



A*B n = 0, 



(61) 



(62) 



separately. The Hermiticity condition, Eq. (|6l|), is very simply satisfied if we 
let 

(63) 



h h 2 



and 



2m' 



(64) 



These relations will be useful in the next section. 
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We now have enough information that we can determine the value of the 
coupling constant, C. Any state in the waveguide must satisfy the condition 
that it be a continuous function of x and y, and that it's slope be a continuous 
function of x and y. Let us consider the state, \^fp). We require that 

(0,y\Q\*p) = (0,v\P\*f>)- (65) 

This, in turn, implies that 

M 

$>j\n(0)&j = X n(0)B n . (66) 
3=1 



If we now compare Eqs. (|64|) and (|66|), we find that the coupling constant is 
given by 

(67) 

In Section (4), we use these results to make contact with scattering theory. 
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Figure 1: The geometry of the two dimensional electron wave guide used in 
our calculations. The rippled waveguide is the region defined with solid lines, 
rectangular waveguide is the region whose upper boundary is given by the dotted 
line. The dotted-dashed line shows the interference between leads and scattering 
region. Here 'a' is the width of the ripple, 'd' is the width of the rectangular 
waveguide, scattering cavity extends from x — to x = L 

Figure 2: Error in cigcnphase, 9, versus energy for an R-Matrix with N = 1000 
and N = 10000 terms kept in the sum and the FEM calculation of 9. The errors 
are calculated in terms of the fractional deviation of the numerically computed 
eigenphase, 8 C , from the exact eigenphase 9 e = 2kL. 

Figure 3: (a) Wigner delay time, t\, for rippled waveguide with parameters, 
a = 10A, d = 100A, L = 200A, for the energies such that only one propagating 
mode exists in the lead. The dashed line is for reaction matrix calculations, the 
solid line is from finite element calculations, (b) Partial Wigner delay times, T\ 
and T2, when there are two propagating modes in the leads. The dashed lines 
are the reaction matrix results and the solid lines are the finite element results. 

Figure 4: Poles of the S-matrix in the lower complex energy plane for an energy 
interval with only one propagating mode in the lead, obtained from Eq. (Jm]). 
The waveguide parameters are a = 10A, d — 100A, L = 200A. 

Figure 5: The effect of evanescent modes for a waveguide with parameters 
a = 10A, d = 100A, L = 500A. (a) The solid line is Wigner delay time, n, 
in the one channel regime, E± < E < 4E\, for an energy interval just before a 
second propagating channel opens in the lead. The dashed line shows T\ when 
no evanescent modes are included in calculations, (b) The amplitudes, ci and c%, 
of the first two evanescent modes in the same energy interval as in (a) . (c) The 
solid line shows the first partial Wigner delay time, t\ just before the opening 
of the third channel. The dashed line is for the case when no evanescent modes 
are included in calculating t\ . (d) The same as (c) for the second partial Wigner 
delay time, T2. 

Figure 6: The absolute value of S-matrix elements, l^nl and 15*121, in two 
modes energy interval, 4£'i < E < 9E\, where two propagating channels are 
allowed. The waveguide parameters are a = 10A, d = 100A, L = 500A. 

Figure 7: Wigner delay times for the energy interval, 4E\ < E < 9Ei, where 
two propagating channels are allowed. The waveguide parameters are a = 10A, 
d = 1001, L = 500A. 
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Figure 8: Poles of |5n | in the lower complex energy plane for the energy interval, 
4Ei < Re(E) < 9Ei. The waveguide parameters are a = 10 A, d — 100 A, L = 
5001. (a) |5n| for < lm(E) < -0.95.Ei. (Only points for which \S n \ < 200 
are shown.) (b) \S n \ for < lm(E) < -0.075.Ei. 



Figure 9: The difference between S-matrix element, \Sn\ with evanescant 
modes, and S-matrix element, \Si±\ without evanescant modes. The waveguide 
parameters are a = 10A, d = 100A, L = 500A. 



Figure 10: The value of |Sii|, when we include the energy dependence of the 
coupling matrices wn p and u>N e (solid line), and its value, |5"° en |, when we fix 
the energy dependence of the coupling matrices Wn p and uiNe to be E = 6.5£i 
(dotted-dashed line). The cavity parameters The waveguide parameters are 
a = 10A, d = 1001, L = 500A. 



Figure 11: The value of \Su\, in the complex energy plane, when we take into 
account the energy dependence of the coupling matrices wn p and WNe (solid 
line), and its value, |6'™ 1 oen |, when we fix the energy dependence of the coupling 
matrices wn p and u>N e to be £ = 6.5Ei (dotted-dashed line). The cavity 
parameters a = 101, d = lOOl, L = 5001. (Only points for which |5n| < 200 
are shown.) 



Figure 12: Histograms of Wigner delay times for different numbers of propa- 
gating channels, M for the case of deterministic scattering ((a)-(c)) and RRMT 
predictions, (d)-(i). (a) Deterministic scattering with M = 2 for a = 10 A, 
d = 1001, L = 500A. (b) Same as (a) for M = 6. (c) Same as (a) for M = 16. 
(d) M = 2 RRMT result with coupling constant, g = 1.8. (e) Same as (d) 
for M = 6. (f) Same as (d) for M = 16. (g) M = 2 RRMT result with a 
unit coupling constant .g = l. (h) Same as (g) for M — 6. (i) Same as (g) for 
M = 16. 
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